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ABSTRACT 

Using a sample of 117053 quasars identified in the PPMXL catalog, the systematic errors of the 
proper motions in the PPMXL catalog are estimated independently as —2.0 and —2.1 mas yr~^ in 
/iaCOS(S and us, respectively. For objects with SDSS r magnitude between 16.0 and 21.0 mag and 
SDSS g — r color between —0.4 and 1.2, there are no obvious magnitude and color dependence in the 
systematic errors of the proper motions in the PPMXL catalog. The random errors of the proper 
motions in the PPMXL catalog change from ~ 4.0 to ~ 8.0 mas yr~^ along with the increase of the 
magnitude of the objects. No color dependence in the random errors of the proper motions in the 
PPMXL catalog is found. There are obvious right ascension a dependence of the systematic errors in 
both components of the proper motions in the PPMXL catalog. For quasars with declination (S > 0, 
the systematic errors of /ia in the PPMXL catalog increase linearly with 5 with a slope of 0.07 mas 
yr~-^ deg~^. 

For comparison, using a subsample of quasars identified both in the PPMXL and SDSS DR7 catalogs, 
the systematic errors of the USNO-SDSS proper motions in the SDSS DR7 catalog are derived as 0.9 
and 0.5 mas yr~^ in /XaCOSiJ and jis, respectively. There are no obvious magnitude, color, and a 
dependence in the systematic errors of the USNO-SDSS proper motions. For quasars with declination 
5 > 0, the systematic errors of in the SDSS DR7 catalog increase linearly with 5 with a slope of 
0.03. The random errors of the USNO-SDSS proper motions change from ^ 2.5 to ^ 6.0 mas yr^^ 
along with the change of the magnitude and color of the objects. Combining the 2MASS or SDSS 
observation data can significantly reduce the systematic errors of the proper motions in the initial 
photographic catalog. 

Subject headings: Stars — Quasars and Active Galactic Nuclei — Data Analysis and Techniques 



1. INTRODUCTION 

Absolute proper motions of stars are very important 
data for the kinematic studies of the Galaxy. For the 
last two decades, great efforts have been made to ac- 
quire high precise and accurate absolute proper motions 
of stars covering the entire sky. The Hipparcos catalog 
is considered as the realization of the International Ce- 
lestial Reference System (ICRS) at optical band. The 
systematic error of the absolute proper motions in the 
Hipparcos catalog with respect to ICRS is estimated to 
be 0.25 mas yr^^. The Hipparcos catalog contains about 
0.12 million star s and is complete down to U = 7.3 — 9.0 
mag (|ESA]|1997n . , 

The Tycho-2 catalog ()Hog et al.ll2"oM ) is the first den- 
sification of the Hipparcos catalog. The proper motions 
in the Tycho-2 catalog are on the system of the Hipparcos 
with a deviation under 0.5 mas yr~^ (Urban et al. 2000). 
The Tycho-2 catalog contains 2.5 million stars and is 
complete down toV ^ 11.5 mag. The ASCC-2.5 catalog 
is based on large, modern, and high-precision catalogs of 
the Hipparcos- Tycho family, including the Tycho-2 cat- 
alog, and provides the most complete all-sky catalog of 
2.5 million stars having un iform h igh-precision astromet- 
ric and photometric data (jKharchenko, .2001.) . 

The UCAC2 catalog presents its proper motions on the 
Hipparcos system with the nominal errors of 1 to 3 mas 
yr~-^ for stars up to U 12 mag and about 4 to 7 mas 
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yr~^ for fainter stars up to y ^ 16 mag. The systematic 
errors of the proper motions in the UC AC2 catalog are 
in th e range from 0.5 to 1.0 mas yr~^ (jZacharias et al.l 
I2004f ). The UCAC2 catalog contains 48 million objects 
with limiting magnitude of i? 16 mag. But UCAC2 
catalog is not a all-sky catalog, most of data in the re- 
gions with declination 5 > +40° do not exist in this 
c atalog. 

iRoser et all (|2008D presented the PPMX catalog which 
including positions and absolute proper motions of ^ 18 
million stars with limi ting magnitude ^ 15 m ag in a red 
band. Very recently, iZacharias et al.l (|2010( ) published 
the all-sky UCAC3 catalog which containing over 100 
million objects with slightly deeper limiting magnitude 
than that of the UCAC2 catalog. 

The absolute proper motions in the above-mentioned 
catalogs can only be used in the kinematic studies of stars 
in the neighborhood of the Sun due to their bright er lim- 
iting magnitudes of stars in those catalogs (e.g. Wu et all 
120091 ). In recent decade, a few of catalogs are avail- 
able for stars with more fainter limiting magnitudes. 
iHamblv et aH (|2001l ) presented the SuperCOSMOS cat- 
alog whose zero-point of the absolute proper motion sys- 
tem is fixed to the extragalactic reference frame. The 
mean error of the proper motions in this catalog changes 
from typically 10 mas yr~^ for objects with i? ~ 18 mag 
to 50 mas yr~^ for objects with i? ~ 21 mag. The zero- 
point error in the proper motions of the SuperCOSMOS 
catalog is less than 1 mas yr^^ for objects with i? > 17 
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mag, and is not larger than 10 mas yr~^ for objects with 
i? < 17 mag. The GSCII catalog presents astrometry, 
photometry, and classification for about 946 million ob- 
jects, but the proper motion data are not released to 
the scientific community due to the significant system- 
atic errors in the derived prop er motions fo r objects in 
the southern hemisphere (L asker et al.ll2008() . 

The USNO-B catalog presents proper motions for 
about one billion objects base d on data obtained from 
scans of 7435 Schmidt plates ()Monet et all 120031 ). But 
the proper motions in this catalog are relative, not ab- 
solute, which limits th eir use in the k inematic studies 
of stars in the Galaxy ()Wu et al.ll2011 f). Combining the 
USNO-B an d Sloan Digit al Sky Survey (SDSS) catalogs, 
IMunn et all ([2003 . [2001 hereafter USNO-SDSS proper 
motion) presented an improved absolute proper-motion 
catalog . The USNO-B positions arc recalibrated us- 
ing the galaxies identified in the SDSS catalog, and the 
proper motions are placed on an absolute reference frame 
and are recomputed using both the USNO-B and SDSS 
positions. The random proper-motion error is only 3 
mas yr~^ for object with r < 18 mag and ~ 5 mas yr~^ 
for object with r < 20 mag. The syst ematic error is typ - 
ically an order of magnitude smaller ([Bond et al.ll2010[ ). 
But the available USNO-SDSS proper motion data are 
li mited by the regions covered by the SDSS observation. 

iFedorov et all ()2010D presented the XPM-1.0 catalog 
which including absolute proper motions of about 280 
millions stars cov ering the entire sk y. The positions in 
the USNO-A2.0 (IMonet et al.l [19981) and 2MASS cata- 
logs were used to calibrate the proper motions of XPM- 
1.0. The zero-point of the absolute proper motions in 
the XPM-1.0 catalog was linked to more than one million 
galaxies. The random error of the proper motions in this 
catalog varies from 3 to 10 mas yr~^ depending on magni- 
tudes. The systematic error of the proper motions in the 
Northern hemisphere is ~ 0.3 mas yr~^, and is ^ 1.0 ma s 
yr~^ in the Southern hemisphere (jFedorov et al.|[2010[ ). 
But, the faint end of XPM-1.0 catalog is limited by the 
limiting magnitude of the 2MASS catalog. 

Very recently, t he PPMXL catalog was released 
(IRoeser et all I2010D. C ombing USNO-B catalog with 
2MASS HSkrutskie et all |2006[) catalog, the new mean 
positions and absolute proper motions of objects listed 
in the USNO-B catalog were derived with least-squares 
adjustment and reduced to the ICRS system using the 
PPMX catalog (Roscr et al. 2008). The PPMXL catalog 
contains ^ 910 millions objects, of which, 410 millions 
objects have photometric data from 2MASS catalog. At 
its bright end, the PPMXL catalog is merged with the 
PPMX catalog. So, PPMXL is the largest collection of 
absolute proper motions on the ICRS system and is com- 
plete from the brightest stars down to U ~ 20 mag. The 
mean errors of the proper motions in PPMXL catalog 
change from 4 mas yr^^ to more than 10 mas yr^^ de- 
pending on the magnitudes. 

As the largest available absolute proper motion cat- 
alog with very fainter limiting magnitude, the PPMXL 
catalog is a very important resource in the kinematic 
studies of the Galaxy. Before using the proper motion 
data in the PPMXL catalog, it is very important to un- 
derstand the systematic a nd random er r ors of the proper 
motions in this catalog. IRoeser et al.l (|2010f) have cor- 
rected the plate-dependent distortions of the proper mo- 



tions du ring the const r uctio n procedure of the PPMXL 
catalog. IRoeser et al.l (|2010[ ) have pointed out that the 
magnitude- and color-dependent systematic errors in the 
PPMXL catalog are difficult to be determined due to no 
independent reference on ICRS at fainter magnitude. On 
the other hand, it is also difficult to distinguish between 
systematic errors and different kinematics for different 
populations in the Galaxy. 

Quasars are very distant extragalactic objects and 
their proper motions can be considered as zero. On the 
other hand, the PSFs of quasars in an image are close 
to those of the stars, the systematic deviation from zero 
within the proper motions of the quasars and their stan- 
dard deviation represent the systematic and random er- 
rors in the proper motions of the stars. In this work, a 
largest available quasars sample is used to analyze the 
systematic and random errors in the proper motions of 
the PPMXL catalog. 

We present the construction procedure of the PPMXL 
catalog in Section 2. The quasar sample identified in the 
PPMXL catalog is presented in Section 3. In section 4, 
the systematic and random errors in the proper motions 
of the PPMXL catalog are derived using the quasar sam- 
ple. For comparison, the systematic and random errors 
in the USNO-SDSS proper motions are derived in Section 
5. Our discussion and conclusions are given in Section 6. 

2. CONSTRUCTION OF THE PPMXL CATALOG 

First, the USNO-B catalog was cross-identified with 
the PPMX catalog. The common stars in the PPMX 
catalog fainter than Tycho-2 limits were used, for each 
field (corresponding to a Schmidt plate), to derive the 
mean deviations in right ascension Aa and declination 
A5 at the field epoch from the USNO-B catalog. This 
step corrected the positions in the USNO-B catalog to 
the ICRS. Then, the field-corrected USNO-B catalog was 
cross-identified with the 2MASS catalog. For each com- 
mon star, its positions at different epochs in the USNO- 
B catalog and at the observational epoch of the 2MASS 
catalog were recalculated with a least-squares method. 
In the least-squares adjustment procedure, individual 
weight for each star at different epoch is attributed. Af- 
ter the least-squares adjustment, a preliminary catalog 
with the new position and proper motion for each com- 
mon star was constructed. 

To link the proper motion system of the preliminary 
catalog to the ICRS, stars in the preliminary catalog with 
Ks between 12 and 13 mag were chosen to represent the 
preliminary system which called PSl, and the PPMX 
catalog was chosen as the optical representation of the 
ICRS. Before linking to the PPMX catalog, the obvious 
plate-dependent distortions of the proper motions of the 
PSl in the south of —20° declination were corrected using 
the proper motions of the UCAC3 catalog. The proper 
motions of the PSl are left unchanged north of —20° 
declination. The UCAC3-corrected PSl system is called 
PS2. Then, using the proper motions in the PPMX cata- 
log, the PS2 was corrected to ICRS, the resulting system 
is called PS3. 

Above-mentioned procedure only considered stars hav- 
ing 2MASS observations. To link the proper motions of 
objects fainter than the 2MASS limits on the ICRS, the 
least-squares adjustment were recalculated giving weight 
zero to the 2MASS observations. For stars in PS3, the 



The systematic and random errors in the PPMXL catalog 



3 



differences in the proper motions with and without in- 
cluding 2MASS observations were derived. Then, these 
differences were used to corrected the proper motions 
of objects without including 2MASS observations to the 
ICRS represented by the PPMX catalog. 

So far, the proper motions of all objects in the PPMXL 
catalog are on the ICRS represented by the PPMX cat- 
alog. At the final stage, PPMXL is merged with PPMX. 
For an object in the PPMX catalog, if there is no object 
matched in the PPMXL or the error of its proper mo- 
tions in the PPMX is smaller than that in the PPMXL, 
the data in the PPMX wiU be added to the PPMXL or 
replace that in the PPMXL. 

3. QUASAR SAMPLE 

The 13th edition of the catalog o f quasars and active 
nuclei ()Veron-Cettv fc VeronI r20ia hereafter VV13) is 
used as the primary source for known quasars. This cat- 
alog contains 133336 quasars which includes all known 
quasars when this catalog was compiled. Two parame- 
ters of quasars: magnitude and coordinate can be used 
to cross-identify the common object between the VV13 
and the PPMXL catalogs. Because the photometric data 
in these two catalogs are photographic magnitude with 
very large uncertainties, only the coordinates of quasars 
in these two catalogs were used. For a quasar listed 
in the VV13 catalog, if only one object is matched in 
the PPMXL catalog within a given search radius, then 
this object is considered as the same one in both two 
catalogs. If the search radius is too small, there is no 
any matched object between these two catalogs. When 
the search radius increases, more objects will be found 
within the given search radius, but no additional data 
can be used to distinguish the true matched one. Us- 
ing a similar method of Souchay ct al. (200§) used in 
their construction of the large quasar astrometric cata- 
log (LQAC), different search radii were used to find the 
one and only matched quasars between the VV13 and 
the PPMXL catalogs. 

Figure [T] shows the number of matched quasars be- 
tween the VV13 and the PPMXL catalogs within differ- 
ent search radii from 0.5 to 30 arcsec with a step of 0.5 
arcsec. It can be seen from Figure [T] that the number 
of matched quasars increases very quickly to its max- 
imum from 0.5 arcsec to 2.5 arcsec, then this number 
decreases slowly along with the increased search radius. 
Using the search radius of 2.5 arcsec, 117053 quasars are 
matched and are taken as the members of our quasar 
sample. For the matched quasars, their data in the PP- 
MXL catalog were extracted using the catalog access 
tool VizieR. Figure [2] shows the spatial distribution of 
the 117053 quasars in our sample. The distribution of 
quasars is very non-uniform, most of quasars are located 
at the limited regions and are spectroscopically con- 
firmed w ith the observations of Sloan Digital Sky Survey 
(SDSS) (lAbazaiian. et al.l [20091: [Schneider et all l2010l) . 
With a search radius of 1 arcsec and using the catalog ac- 
cess tool VizieR, the coordinates of quasars listed in the 
PPMXL catalog were used to cross-identifi ed with the 
seven th data release (DR7) of SDSS lAbazaiian. et al.l 
120091) ■ 97150 quasars in our sample were identified in the 
SDSS DR7 catalog. 

Figures [3] and [4] show the magnitude and color distri- 
bution of quasars in our sample, respectively. The B 



and R magnitudes used in Figures [3] and |3| are photo- 
graphic magnitudes derived from the PPMXL catalog. 
The shadow regions in the left panels of Figures [3] and 
[4| represent the distribution of the quasars whose SDSS 
data are available. The SDSS r magnitude and g — r color 
distributions for the quasars in our sample are plotted in 
the right panels of Figures [3] and (H respectively. Fig- 
ure [3] shows that most of quasars in our sample have the 
photographic B magnitudes between 17.0 and 21.5 mag 
and have SDSS r magnitudes between 17.0 and 20.5 mag. 
Figure [3| indicates that our present quasar sample repre- 
sents the objects with fainter magnitudes in the PPMXL 
catalog. In our quasar sample, only 13520 quasars have 
the observation data from the 2MASS catalog. Thus, 
most of the proper motions of quasars in our sample 
were derived only from the photographic data. Figure 
[4| shows that most of quasars in our sample are within 
the range oi B — R between —1.0 and 2.0 and of SDSS 
g — r between —0.2 and 0.8. Figure [3| also indicates that 
our present quasar sample only represents the objects in 
the PPMXL catalog within the limited color range. 

4. THE SYSTEMATIC AND RANDOM ERRORS OF THE 
PROPER MOTIONS IN THE PPMXL CATALOG 

We have assumed that a quasar's distance is very far 
that its proper motion cannot be detected by present 
astrometric observation in the optical band and can be 
considered as zero. Thus, for a sample of quasars, any 
deviation from zero of the mean of the proper motions 
of the sample can be considered as the systematic error 
of the proper motions, and the dispersion of the proper 
motions from the mean can be considered as the random 
error of the proper motions of the sample. Figure [5| shows 
the histograms of the proper motions of the quasars in 
our sample in right ascension cos 5 and in declination 
/i5 with a bin of 1.0 mas yr~^. Figure [S] only shows the 
distributions of the proper motions for quasars whose 
proper motions with the absolute value less than 35 mas 
yr""'^ in each component. A Gaussian function was used 
to fit the distribution of the proper motions in each com- 
ponent of all quasars in our sample and was plotted as 
red solid line in each panel of Figure O The region under 
the red solid line in each panel includes proper motion 
data within the 3(t of the mean. The best-fitting pa- 
rameters of the Gaussian function: the mean and the 
dispersion of the proper motions in each component are 
labeled in each panel of Figure [S] Figure [5] indicates that 
the mean or the systematic error of the proper motions 
of quasars in the PPMXL catalog is —2.0 and —2.1 mas 
yr~^ in /Iq cos (5 and /x^, respectively. The random error 
of the proper motions of quasars in the PPMXL catalog 
is 5.3 and 4.8 mas yr~^ in ^aC0s5 and us, respectively. 

4.1. The magnitude dependence of the systematic and 
random errors of the proper motions 

Figure [Bj shows the magnitude dependence of the 
proper motions of quasars in our sample in each compo- 
nent. The top two panels are based on the photographic 
B magnitude, and the bottom two panels are based on 
the SDSS r magnitude. The blue dashed line in each 
panel shows the proper motions with value of zero. In 
order to derive the dependence of the proper motions 
in each component on the magnitude, quasars are di- 
vided into different magnitude groups with a bin of 0.5 
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mag. For photographic B magnitude, quasars with B 
magnitude in the range between 16.0 and 22.0 mag were 
used. For SDSS r magnitude, quasars with r magnitude 
in the range between 16.0 and 21.0 mag were used. For 
each magnitude group of quasars, a Gaussian function 
was used to fit the proper motions in each component in 
that group, the mean and dispersion of the best-fitting 
Gaussian function corresponding to the systematic and 
random errors of proper motions of quasars in that group 
are showed as red open circle and error bar in each panel 
of Figure [6l 

The left panels of Figure |6] indicate that the system- 
atic errors of /Xq, cos S have the similar dependence on 
the photographic B magnitude and on the SDSS r mag- 
nitude. The absolute value of the systematic deviation 
of fia cos d increases from the brightest end of the quasar 
magnitude group to its maximum, then decreases pro- 
gressively to the faintest end of the quasar magnitude 
group. The rms of the systematic errors of Ha cos 6 for 
different quasar magnitude groups is 0.7 mas yr~^ for 
photographic B magnitude and 0.5 mas yr~^ for SDSS r 
magnitude, respectively. But this rms is smaller than the 
absolute systematic deviation of 2.0 mas yr^^ in /Iq, cos S 
for the whole quasar sample by a factor of 3 - 4. 

Except the last two points in the faint photographic 
B magnitude bins, the right panels of Figure H] indicate 
that the systematic errors of fis have the similar depen- 
dence on the photographic B magnitude and on the SDSS 
r magnitude. The absolute value of the systematic de- 
viation of Us decreases progressively from the brightest 
end of the quasar magnitude group to the fainter end 
of the quasar magnitude group. The rms of the system- 
atic errors of /ia for different quasar magnitude groups 
is 0.5 mas yr~^ for photographic B magnitude and 0.4 
mas yr~^ for SDSS r magnitude, respectively. This rms 
is smaller than the absolute systematic deviation of 2.1 
mas yr~^ in for the whole quasar sample by a factor 
of 4. 

The red error bar in each panel of Figure |6] shows the 
random error of the proper motions in each component 
in each quasar magnitude group. All of the four panels 
of Figure |6] indicate that the random errors of proper 
motions in each component increase along with the mag- 
nitude of quasars. The magnitude dependence of the 
random error of proper motions in each component can 
be described by a function: cr^ = a + bx exp'^™~^\ where 
cr^ is the random error of proper motions in one com- 
ponent, a, b, and c are the unknown parameters to be 
fitted, and m is the B ov r magnitude. In each panel of 
Figure El the yellow dash dot line shows the best-fitting 
function. Table [T] lists the best-fitting parameters for dif- 
ferent proper motions components and magnitude types. 
In general, the random error of proper motions in one 
component changes from ~ 4.0 mas yr~^ to ~ 8.0 mas 
yr^^ along with the magnitude changing from 16.0 to 
21.0 mag. 

4.2. The color dependence of the systematic and 
random errors of the proper motions 

Figure [7] shows the color dependence of the proper mo- 
tions of quasars in our sample in each component. The 
top two panels are based on the photographic B — R 
color, and the bottom two panels are based on the SDSS 
g — r color. The blue dashed line in each panel shows the 



proper motions with value of zero. In order to derive the 
dependence of the proper motions in each component on 
the color, quasars are divided into different color groups 
with a bin of 0.5 mag ioi B — R and a bin of 0.2 mag 
for g — r, respectively. For photographic B — R color, 
quasars with B — R'm the range between —1.25 and 2.75 
were used. For SDSS g — r color, quasars with g — r in the 
range between —0.4 and 1.2 were used. For each color 
group of quasars, a Gaussian function was used to fit the 
proper motions in each component in that group, the 
mean and dispersion of the best-fitting Gaussian func- 
tion corresponding to the systematic and random errors 
of proper motions of quasars in that group are showed as 
red open circle and error bar in each panel of Figure [T] 

Figure [7] indicate that there is no obvious color depen- 
dence of the systematic errors of fia cos d and fis specially 
in the g — r panel. The rms of the systematic errors in 
^laCOsS and for different quasar color groups is 0.8 
and 0.3 mas yr~^ for the photographic B — R color, re- 
spectively; and 0.2 and 0.4 mas yr~^ for the SDSS g — r 
color, respectively. In general, the dispersion of system 
errors of proper motions in one component for quasars in 
different color group is smaller than the absolute system- 
atic deviation of the proper motions in that component 
for all quasars in the sample by a factor of 2 - 4. 

Figure [7] indicates that there is no obvious color de- 
pendence of the random errors of proper motions in each 
component both for B — R and g — r. But for quasars with 
very bluer or redder color, the random errors of proper 
motions of them are bigger than those for quasars with 
color in the range of —0.2 < g — r < 0.7 . The rms of 
the random errors of the proper motions for quasars in 
different color groups is 0.8 mas yr^^ in fia cos (5 and 0.5 
mas yr"^ in fig, respectively. This rms is much smaller 
than the random error of the proper motions in each 
component for all quasars in our sample by a factor of 5 
- 10. 

4.3. The position dependence of the systematic and 
random errors of the proper motions 

Figure [8] shows the proper motions of quasars in our 
sample in each component as a function of right ascension 
a and declination S . The blue dashed line in each panel 
shows the proper motions with value of zero. In order 
to derive the dependence of the proper motions in each 
component on the position, quasars were divided into dif- 
ferent position groups with a bin of 35° in a and a bin of 
15° in S, respectively. For each position group of quasars, 
a Gaussian function was used to fit the proper motions 
in each component in that group, the mean and disper- 
sion of the best-fitting Gaussian function corresponding 
to the systematic and random errors of proper motions 
of quasars in that group are showed as red open circle 
and error bar in each panel of Figure [H 

The top two panels of Figure [5] shows the distribu- 
tion of fia cos i5 and for quasars in our sample along 
with a. These panels of Figure [8] indicate that, both 
for fia cos S and fis, there are systematic dependence of 
the systematic errors of the proper motions on a. A 
function oi Ji — a + bx sin(a — c) was used to fit the 
a dependence of the proper motions in each compo- 
nent, where JI is the systematic error of the proper mo- 
tions in each component for quasars in each a bin, a, 
b, and c are the unknown parameters to be fitted, and 
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a is in degree. For /x^ cos S, the best-fitting result of 
the function is fiaCOsS = —0.59 + 2.95sin(a + 78.48). 
For iJs, the best- fitting result of the function is Jig — 
-2.30 - 0.99 sin(a -I- 61.02). The best-fitting function for 
the proper motions in each component is plotted as yel- 
low dash dot line in each top panel of Figure El The 
top two panels of Figure |8] also indicate that there is no 
systematic dependence of the random error of the proper 
motions in each component on a. The dispersion of the 
random errors in /ic cos S and fj,s for quasars in different 
a bins is 1.0 and 0.8 mas yr~^, respectively. This disper- 
sion is smaller than the random error of proper motions 
for all quasars in our sample by a factor of 4 - 6. 

The bottom two panels of Figure |8] shows the distri- 
bution of fia cos S and fig along with S. The bottom-left 
panel of Figure [5] indicates that, for quasars with 6 > 0, 
the systematic errors of /Iq cos S in the six S bins are 
smaller than zero. On the other hand, for quasars with 
S < 0, the systematic errors of /Xq cos d in the four d bins 
are bigger than zero. For quasars with 5 > 0, there is 
no S dependence of the systematic and random errors in 
/Zq cos(S, the dispersion of the systematic and random er- 
rors of iia cos (5 for quasars in different S bins is 0.5 and 
0.7 mas yr~^, respectively. For quasars with ^ < 0, the 
maximum systematic error of cos (5 = 4.6 mas yr~^ is 
in the nearby of (5 = —30°. 

The bottom-right panel of Figure [8] indicates that, for 
quasars with 6 < 0, there is no S dependence of the 
systematic and random errors of /i^, but the quasars in 
the nearby of d — —30° have the maximum systematic 
error of Jig — —5.3 mas yr^^. For quasars with (5 > 0, 
the bottom-right panel of Figure [8] indicates that there 
is S dependence of the systematic error of fis- A line 
of Jig = —3.65 -I- 0.07(5 can best fit the S dependence of 
^s, where ]Is is the systematic error of fis for quasars in 
different S bins, 6 is in degree. This line is plotted as 
yellow dash dot line in the bottom-right panel of Figure 
m For quasars with S > 0, there is no 5 dependence of the 
random errors of fig , the dispersion of the random errors 
of for quasars in different S bins is 0.3 mas yr~^. 

4.4. The importance of the 2MASS observation data 

In the previous sections, we have pointed out that there 
is only 13520 or 12% quasars in our sample with the 
2MASS observation data. Thus, the proper motions for 
most of quasars in our sample were derived only from 
the photographic data. But it is important to compare 
the derived proper motions with and without the 2MASS 
data to understand the importance of the 2MASS data 
in derivation of the proper motions in the PPMXL cata- 
log. A appropriate quasar sample was chosen before this 
comparison. The top two panels of Figure [9] show the 
photographic B magnitude distributions of quasars with 
(left panel) and without (right panel) the 2MASS data. 
These two panels indicate that the B magnitude distri- 
butions for these two groups of quasars are very different. 
Most of quasars with 2MASS data are bright and the me- 
dian of the B magnitude distribution is ^ 18.0 mag. On 
the other hand, most of quasars without 2MASS data are 
faint and the median of the B magnitude distribution is 
~ 20.0 mag. Furthermore, for quasars with B > 19.0, the 
number of quasars without 2MASS data is much bigger 
than that of quasars with 2MASS data, and all quasars 
without 2MASS data whose B magnitudes are bigger 



than 17.0 mag. Thus, quasars with B magnitudes in 
the range between 17.3 and 18.3 mag were chosen as the 
comparison sample. The number of quasars with and 
without the 2MASS data in the comparison sample is 
6013 and 5878, respectively. 

The middle and bottom panels of Figure [3] show the 
histograms of the proper motions in each component for 
quasars with (left panels) and without (right panels) the 
2MASS data in the comparison sample , respectively. 
Similar to Figure [SJ A Gaussian function was used to fit 
the proper motions in each component and was plotted 
as red solid line in each panel of Figure [HI These four 
panels in Figure [HI indicate that there is no obvious dif- 
ference in the random errors of the proper motions in 
both components for quasars with and without 2MASS 
data. But the systematic error of fXa cos S for quasars 
with 2MASS data is much smaller than that for quasars 
without 2MASS data. The systematic error of /i^ for 
quasars with 2MASS data is also a little smaller than 
that for quasars without 2MASS data. Figure [9] indi- 
cates that the 2MASS data are very important in the 
derivation of the proper motions in the PPMXL catalog. 
They can significantly reduce the systematic errors of the 
derived proper motions in the PPMXL catalog. 

5. COMPARISON WITH THE USNO-SDSS PROPER 
MOTION 

The USNO-SDSS proper motion uses the USNO-B cat- 
alog as the early echo data, but links its zero-point to the 
galaxies. The SDSS photometric data as the later echo 
data used in deriving the USNO-SDSS proper motion are 
much deeper than the 2MASS photometry used by the 
PPMXL catalog. 83% quasars in the PPMXL catalog 
are also found in the SDSS DR7. Thus, it is important 
to compare the proper motions of quasars in the PPMXL 
catalog to that in the SDSS catalog to see whether the 
new SDSS photometric data used in the derivation of 
the USNO-SDSS proper motion can improve the derived 
proper motions for the given quasar sample. 

Similar to Figure \E\ Figure [10] shows the histograms 
of the USNO-SDSS proper motions of quasars cross- 
identified in the PPMXL and SDSS DR7 catalogs with 
a bin of 1.0 mas yr^^. A Gaussian function was used 
to fit the USNO-SDSS proper motions in each compo- 
nent and was plotted as red solid line in each panel of 
Figure [TOl Figure [TOj indicates that the systematic er- 
ror of the USNO-SDSS proper motions of quasars in the 
SDSS DR7 catalog is 0.9 and 0.5 mas yr^^ in /i^ cos 6 and 
fis, respectively. The random error of the USNO-SDSS 
proper motions of quasars in the SDSS DR7 catalog is 3.3 
and 3.3 mas yr~^ in fiaCosS and fig, respectively. The 
derived systematic and random errors are much smaller 
than those of the proper motions of quasars in the PP- 
MXL catalog, which indicates the improvement in the 
derivation of proper motions including the SDSS astro- 
metric and photometric data. 

Similar to Figure [6l Figure [Til shows the SDSS r mag- 
nitude dependence of the USNO-SDSS proper motions 
in each component of quasars cross-identified in the PP- 
MXL and SDSS DR7 catalogs. Quasars with SDSS r 
magnitude in the range between 16.0 and 21.0 mag were 
divided into different magnitude groups with a bin of 
0.5 mag. For each magnitude group of quasars, a Gaus- 
sian function was used to fit the proper motions in each 
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component in that group. Figure [TT] indicates that there 
is no obvious magnitude dependence of the systematic 
errors of the USNO-SDSS proper motions in the both 
two components. The rms of the systematic errors of 
the USNO-SDSS proper motions for quasars in different 
magnitude groups is 0.3 and 0.4 mas yr~^ in /ZqCOsJ 
and fis, respectively. This rms is a httle smaller than 
that for quasars in different magnitude groups in the 
PPMXL catalog. Figure [TT] also indicates that the ran- 
dom errors of the USNO-SDSS proper motions in each 
component increase along with the SDSS r magnitude of 
quasars. Similar to Figure El the yellow dash dot line in 
each panel of Figure [TT] shows the best-fitting function 
which describing the SDSS r magnitude dependence of 
the random error of the USNO-SDSS proper motions in 
each component. The best-fitting function for the ran- 
dom errors in /Xa cos 5 is cr^^ cos 5 — 2.3 -I- 3.3 exp^''^^"-^) 
mas yr~^, and the best-fitting function for the random 
errors in fig is (t^^ — 2.3 -f- 3.9 exp(''~^°-*) mas yr^^. In 
general, the random error of the USNO-SDSS proper mo- 
tions in one component changes from ~ 2.5 mas yr~^ to 
~ 6 mas yr~^ along with the SDSS r magnitude changing 
from 16.0 to 21.0 mag. 

Similar to Figure [T] Figure [12] shows the SDSS g - r 
color dependence of the USNO-SDSS proper motions in 
each component of quasars cross-identified in the PP- 
MXL and SDSS DR7 catalogs. Quasars with SDSS g-r 
color in the range between— 0.4 and 1.4 were divided into 
different color groups with a bin of 0.2 mag. For each 
color group of quasars, a Gaussian function was used to 
fit the proper motions in each component in that group. 
Figure [T^ indicates that there is no obvious color depen- 
dence of the systematic errors of the USNO-SDSS proper 
motions in the both two components. The rms of the 
systematic errors of the USNO-SDSS proper motions for 
different quasar color groups is 0.1 and 0.3 mas yr^^ in 
cos (5 and fis, respectively. This rms is a little smaller 
than that for quasars in different color groups in the PP- 
MXL catalog. Figure [T^ also indicates that the random 
errors of the USNO-SDSS proper motions in both two 
components increase along with the SDSS g — r colors of 
quasars. The color dependence of the random errors of 
the USNO-SDSS proper motions in each component was 
fitted and showed as yellow dash dot line in each panel 
of Figure [T^] using the similar function used in Figures 
[1] and [n] (T^ = a -f 6 X exp^™-'^), here m is the SDSS 
g — r color. The best-fitting function for the random 



(m-1.5) 



mas 



errors in /i^ cos 5 is <7^„cos5 = 2.4 -I- 3.9 exp 
yr~^, and the best-fitting function for the random errors 
in is (Jfj.s = 2.3 -I- 3.7exp("'~^-^) mas yr~^. In general, 
the random error of the USNO-SDSS proper motions in 
one component changes from ~ 2.5 mas yr^^ to 6 mas 
yr""'^ along with the SDSS g — r color changing from —0.5 
to 1.5. 

Similar to Figure [S] Figure [H shows the USNO-SDSS 
proper motions in each component of quasars cross- 
identified in the PPMXL and SDSS DR7 catalogs as a 
function of a and S. The top two panels of Figure [13] 
indicate that there is no a dependence of the systematic 
and random errors of the USNO-SDSS proper motions 
of quasars in the both two components. The rms of the 
systematic errors of the USNO-SDSS proper motions for 
quasars in different a groups is 0.2 and 0.4 mas yr^^ in 



/^Q. cos 6 and fis , respectively. The rms of the random 
errors of the USNO-SDSS proper motions for quasars in 
different a groups is 0.4 and 0.4 mas yr~^ in /ZqC0S(5 
and lis, respectively. The bottom- left panel of Figure 
1131 indicates that there is no obvious S dependence of 
the systematic and random errors of the USNO-SDSS 
proper motions in cos S. The rms of the systematic 
and random errors of the USNO-SDSS proper motions 
in cos S for quasars in different i5 groups is 0.2 and 0.4 
mas yr~^, respectively. The bottom-right panel of Fig- 
ure [13] indicates that there is no obvious S dependence 
of the random errors of the USNO-SDSS proper motions 
in ^s- The rms of the random errors of the USNO-SDSS 
proper motions in /j^g for quasars in different 6 groups is 
0.4 mas yr~^. Similar to the bottom-right panel of Fig- 
ure [8] the bottom- right panel of Figure [T3l indicates the S 
dependence of the systematic errors of the USNO-SDSS 
proper motions in ^s- For quasars with (5 > 0, a line of 
JIs = —0.18-1-0.03(5 was used to best fit the i5 dependence 
of the systematic errors of fxs and showed as yellow dash 
dot line in the bottom- right panel of Figure [13] The co- 
efficient 0.03 is half of that derived for the S dependence 
of the systematic errors of us for quasars in the PPMXL 
catalog. 

6. SUMMARIES 

As the very distant extragalactic source, the proper 
motion of a quasar has been assumed to be zero. Proper 
motions of quasars with very high accuracy are needed 
to check this assumption. In order to detect the sec- 
ular aberration drift of the extragalactic radio source 
proper motions caused by the rotation of th e Solar Sys- 
tem b arycenter around the Galactic center, iTitov et al.l 
()2011t ) derived the absolute proper motions of 555 ex- 
tragalactic radio sources based on the observation by 
the very long baseline interferometry (VLBI) between 
1990 and 2010. For the most observed radio sources, 
the i nflated position e rrors are 40 /xas in each coordi- 
nate ()Titov et al.ll2011l). The mos t of extragalactic radio 
sources used bv lTitov eTaLl ([20T1II are quasars , thus, the 
proper motions derived bv iTitov et al.l (|201lD are those 
with the highest accuracy for quasars at present. 

Figure [T^ shows the h i stogra ms of the proper motions 
derived bv lTitov et all (|20ll in fj,a cos 6 and fis with 
a bin of 5 uas yr ~^. Not all of 555 sources used by 
ITitov et al.l ([201I are showed in Figure I14[ and sources 
with the absolute value of proper motion in each com- 
ponent bigger than 1 mas yr"-'^ are not showed in Figure 
[T4l The number of sources not showed in Figure [14] is 13 
and 7 in the fia cos 5 and the panel, respectively. Sim- 
ilar to Figure [5] a Gaussian function was used to fit the 
proper motio ns in each compone nt for all 555 sources in 
the sample of ITitov et al.l ()201lD and was plotted as red 
solid line in each panel of Figure [14] Figure [141 indicates 
that t he mean of proper motions derived bv lTitov et al.l 
(|20Tl[) for 555 sources is 0.0 and —0.9 /ias yr ^ in /i^ cos S 
and /i5 , respectively. The rm s of the derived proper mo- 
tions by ITitov et al.l (poTTl ) for 555 sources is 27.9 and 
35.4 /ias yr~^ in /i^cos^ and fig, respectively. It should 
be noted that the unit in Figure [14] is /las yr~^ not mas 
yr~^ used in Fi gures 15] and ITUl 70% sources in the sample 
of lTitov et all |20T1 have proper motions with values in 
the range within the 3cr of the mean in each component. 

Figure [14] indicates that the mean of proper motions 
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with high accuracy for quasa rs is cfose to zero. Based on 
a sample of 62 radio sources, I Moor et all ()201lD pointed 
out that the proper motions of the radio sources are 
caused by the changes in their brightness distribution 
structure. The systematic motion of quasars can be de- 
tected at present is the sec ular aberration dr ift with a 
value of 6.4 ± 1.5 ^ias yr"! (jTitov et al.ll20lH) . In gen- 
eral, the mean of proper motions of quasars can be con- 
sidered as zero with present observation accuracy. The 
systematic motion of 6.4 ^as yr~^ of quasars cannot be 
detected by the present ground-based observation in the 
optical band. 

Based on the assumption that the mean of proper mo- 
tions for a sample of quasars should be zero, any devia- 
tion from zero of the mean of the proper motions of the 
sample can be considered as the systematic error of the 
proper motions, and the dispersion of the proper mo- 
tions from the mean can be considered as the random 
error of the proper motions of the sample. A sample of 
117053 quasars identified in the PPMXL catalog were 
used to derive the systematic and random errors in the 
proper motions of the PPMXL catalog. The systematic 
error of the proper motions for the whole quasar sample 
is derived as —2.0 and —2.1 mas yr"-'^ in Ha cos 6 and ^s, 
respectively. The random error of the proper motions 
for the whole quasar sample is 5.3 and 4.8 mas yr^^ in 
Ha cos S and , respectively. 

No obvious magnitude dependence of the systematic 
errors of the proper motions in both components are 
found. The rms of the systematic error of the proper 
motions in each component for quasars with different 
magnitudes is ~ 0.5 mas yr~^. There is obvious mag- 
nitude dependence of the random errors of the proper 
motions in both components. The random error of the 
proper motions in each component can change from ~ 4.0 
to ~ 8.0 mas yr^^ along with the increase of the mag- 
nitude of quasars. No obvious color dependence of the 
systematic and random errors of the proper motions in 
both components are found. The rms of the system- 
atic errors of the proper motions in each component for 
quasars with different colors is ^ 0.5 mas yr~^. The rms 
of the random error of the proper motions in each com- 
ponent for quasars with different colors is ^ 0.8 mas yr^^ 
and ~ 0.5 mas yr~^ for photographic B — R and SDSS 
g — r, respectively. 

The systematic errors of the proper motions in both 
components Ji are correlated with a of quasars: Jl oc 
sin a. No obvious a dependence of the random errors of 
the proper motions in both components are found. The 
rms of the random error of the proper motions in each 
component for quasars with different a is ~ 1.0 mas yr~^. 
For quasars with 5 < 0, their random errors of the proper 
motions in both components are bigger than those for 
quasars with (5 > 0; and their systematic errors of the 
proper motions in both components are very different 
from those for quasars with 5 > Q. Especially, quasars 
in the neighborhood oi 5 = —30° have the maximum 
of the systematic error of the proper motions in both 
components ~ 5 mas yr"'^. For quasars with (5 > 0, 
there is no obvious 6 dependence of the systematic and 
random errors of cos 5. The rms of the systematic and 
random errors of cos 5 for quasars with 5 > is 0.5 
and ^ 0.7 mas yr~^, respectively. But for quasars with 
5 > 0, the systematic errors of us increase with the 5 and 



can be best described as /J^ oc 0.07i5. There is no obvious 
5 dependence of the random errors of us for quasars with 
(5 > 0. The rms of the random errors of /i^ for quasars 
with (5 > is ^ 0.3 mas yr~^. 

The 2MASS observation data are important in the 
derivation of the proper motions in the PPMXL catalog. 
The systematic error of cos 6 for quasars with 2MASS 
data is about half of that for quasars without 2MASS 
data. For comparison, the systematic and random errors 
of the USNO-SDSS proper motions for a quasar sample 
identified in the PPMXL and SDSS DR7 catalogs are de- 
rived. The systematic error of the USNO-SDSS proper 
motions is 0.9 and 0.5 mas yr~^ in ^qC0S(5 and ^s, re- 
spectively. The random error of the USNO-SDSS proper 
motions is 3.3 and 3.3 mas yr~^ in fiaCOsS and /x^, re- 
spectively. 

No obvious magnitude dependence of the systematic 
errors of the USNO-SDSS proper motions in both com- 
ponents arc found. The rms of the systematic errors of 
the USNO-SDSS proper motions in each component for 
quasars with different magnitudes is ^ 0.4 mas yr^^. 
The random error of the USNO-SDSS proper motions in 
each component increases from 2.5 mas yr^^ to 6.0 
mas yr^^ along with the increase of the SDSS r mag- 
nitude of quasars. No obvious color dependence of the 
systematic errors of the USNO-SDSS proper motions in 
both components is found. The rms of the systematic 
errors of the USNO-SDSS proper motions in each com- 
ponent for quasars with different g — r color is ~ 0.3 mas 
yr~^. But, the random errors of the USNO-SDSS proper 
motions in each component can change from '^2.5 mas 
yr"-'^ to ~ 6.0 mas yr~^ along with the SDSS g — r color 
of between —0.5 and 1.5. 

There is no a dependence of the systematic and ran- 
dom errors of the USNO-SDSS proper motions in both 
two components. The rms of the systematic and random 
errors of the USNO-SDSS proper motions in each com- 
ponent for quasars with different a is ^ 0.4 and ~ 0.4 
mas yr~^, respectively. There is no 5 dependence of 
the random errors of the USNO-SDSS proper motions 
in both two components. The rms of the random errors 
of the USNO-SDSS proper motions in each component 
for quasars with different 5 is ~ 0.4 mas yr^^. No obvi- 
ous 5 dependence of the systematic errors of the USNO- 
SDSS proper motions in cos 5 is found. The rms of the 
systematic errors of the USNO-SDSS proper motions in 
/iaCOSfJ for quasars with different (5 is ^ 0.2 mas yr^^. 
The systematic errors of the USNO-SDSS proper mo- 
tions in jjLs increase with the 5 and can be best described 
as oc 0.035. 

In general, comparing with the proper motions in the 
PPMXL catalog, using the SDSS photometric and as- 
trometric data can significantly reduce the systematic 
and random errors of the derived proper motions. Thus, 
combining the photographic data and the present obser- 
vation data obtained from the ground-based telescopes 
with CCD, the proper motions with systematic error less 
than 1 mas yr~^ and random error less than 5 mas yr~^ 
for objects with r < 20 mag can be derived. 
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Fig. 1. — Number of the uniquely cross-identified quasars between the W13 and the PPMXL catalogs within the given search radius. A 
search radius of 2.5 arcsec was adopted finally. 




Fig. 2. — Spatial distribution of 117053 quasars in the PPMXL catalogs in a J2000 equatorial reference frame with a Aitoff projection. 



The systematic and random errors in the PPMXL catalog 



9 




5.0x10 



14 16 18 20 22 24 
B (mog) 




14 16 18 20 
r (mog) 



22 



Fig. 3. — Magnitude distributions of quasars in the PPMXL catalog. Left panel: photographic B magnitude distribution. The shadow 
area is the distribution of cross-identified quasars between the PPMXL and the SDSS DR7 catalogs. Right panel: SDSS r magnitude 
distribution. 
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Fig. 4. — Color distributions of quasars in the PPMXL catalog. Left panel: photographic B — R color distribution. The shadow area is 
the distribution of cross-identified quasars between the PPMXL and the SDSS DR7 catalogs. Right panel: SDSS g — r color distribution. 
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Fig. 5. — Histograms showing the proper motion distributions of quasars in the PPMXL catalog. The red solid line in each panel is the 
best-fitting Gaussian function. The best-fitting parameters of the Gaussian function: the mean and the dispersion are also labeled in each 
panel. Only quasars with the proper motions in each component with the absolute value less than 35 mas yr~^ are showed in this figure, 
but the proper motions for all quasars in our sample in each component are used to fit the Gaussian function. The region under the red 
line in each panel includes proper motion data within the 3cr from the derived mean. 



TABLE 1 
The best-fitting parameters of 
FUNCTION: (T^ = a + fe X exp('"~'^) 
DESCRIBING THE MAGNITUDE DEPENDENCE 
OF THE RANDOM ERRORS OF PROPER 
MOTIONS IN EACH COMPONENT DERIVED 
FROM THE QUASARS IN THE PPMXL 
CATALOG, a, b, AND C ARE UNKNOWN 
PARAMETERS TO BE FITTED. 
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type of (T^ type of m a 



cos 5 
cos 5 

IJ-s 



B 



B 



4.4 3.4 21.0 

3.8 3.3 20.5 

4.2 2.1 21.0 

3.7 2.6 20.5 
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Fig. 6. — Magnitude dependence of the proper motions of quasars in the PPMXL catalog. Top panels: photographic B magnitude. 
Bottom panels: SDSS r magnitude. The blue dashed line in each panel shows the proper motions with value of zero. The yellow dash 
dot line in each panel is the test-fitting function which describes the magnitude dependence of the random errors of proper motions in the 
PPMXL catalog: cr^ = a + b X exp'"'"'^) , where is the random error of proper motions in one component, a, b, and c are the unknown 
parameters to be fitted, and m is the B or r magnitude. Not all of the quasars in our sample are showed, but the proper motions in each 
component for all quasars in each magnitude bin are used to derive the systematic and random errors in that magnitude bin, which are 
showed as red open circles and error bars in each panel. 
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Fig. 7. — Color dependence of the proper motions of quasars in the PPMXL catalog. Top panels: photographic B — R color. Bottom 
panels: SDSS g — r color. The blue dashed line in each panel shows the proper motions with value of zero. Not all of the quasars in 
our sample are showed, but the proper motions in each component for all quasars in each color bin are used to derive the systematic and 
random errors in that color bin, which are showed as red open circles and error bars in each panel. 
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Fig. 8. — Proper motions of quasars in the PPMXL catalog as a function of a and 5. The blue dashed line in each panel shows the proper 
motions with value of zero. Not all of the quasars in our sample are showed, but the proper motions in each component for all quasars 
in each position bin are used to derive the systematic and random errors in that position bin, which are showed as red open circles and 
error bars in each panel. The yellow dash dot lines in the top panels show the best- fitting function: Ji = a + b X sin(o — c), where Ji is the 
systematic error of proper motions in each component of quasars in each a bin, a, b, and c are the unknown parameters to be fitted, and a 
is in degree. The yellow dash dot line in the bottom-right panel is the best-fitting line for quasars with 5 > 0: Jig = —3.65 + 0.07<5, where 
/JJ is the systematic error of /ig of quasars in each S bin, 5 is in degree. 



14 



Wu et al. 



4000 



3000 



a: 2000 



1000 




600 
500 
400 
^ 300 
200 

100 




16 


18 20 


22 




B (mog) 






n /XaCOS(5(j = 


-1.4 




a: 


4.4 \ 









2.5x1 0"! ^ 

2.0x10'' ■ 
1.5X10'' ; 
1.0x10* ■ 

5.0x10^ - 

0[ i_ 



16 



18 20 
B (mog) 



22 




-30 -20 -10 10 20 30 
fXgCOsd (mos yr"') 



-30 -20 -10 10 20 30 
^„cos<5 (mos yr"') 




-20 -10 10 20 
M« (mos yr"') 




-10 10 20 
Ma (mos yr"') 



Fig. 9. — Histograms of the photographic B magnitudes and the proper motions of quasars with and without 2MASS observation data 
in the PPMXL catalog. Left panels are for quasars with 2MASS data, and right panels are for quasars without 2MASS data. Only quasars 
with B in the range between 17.3 and 18.3 mag are used in the middle and bottom panels showing the proper motions distributions. 
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Fig. 10. — Similar to Figure [S] histograms of the USNO-SDSS proper motions for quasars identified in the PPMXL and SDSS DR7 
catalogs. 
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Fig. 11. — Similar to Figure[6l SDSS r magnitude dependence of the USNO-SDSS proper motions for quasars identified in the PPMXL 
and SDSS DR7 catalogs. 
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Fig. 12. — Similar to FigureFfl SDSS g — r color dependence of the USNO-SDSS proper motions for quasars identified in the PPMXL and 
SDSS DR7 catalogs. The yellow dash dot line in each panel is the best-fitting function which describing the color dependence of the random 
errors of the USNO-SDSS proper motions: = a + b X exp'™"'^', where cr^ is the random error of proper motions in one component, a, 
b, and c are the unknown parameters to be fitted, and m is the SDSS g — r color. 
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Fig. 13.— Similar to Figure[8l but for the USNO-SDSS proper motions for quasars identified in the PPMXL and SDSS DR7 catalogs. 
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Fig. 14. — Si milar to Figure [5] histograms of the proper motions of extragalactic sources derived with the observation by VLBI between 
1990 and 2010 ( Titov et al.ll201ir) . Not all 555 sources are showed in this figure, the number of outliers is 13 and 7 in the fi^ cos (5 and the 
fis panel, respectively. 



